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III. REVIEW OF BASIC Quantum Mechanics: Two- Level 

Quantum Systems: 

The literature of quantum optics and laser spectroscopy abounds with discussions of the 
two-level (two-state) system. This emphasis comes about because the interaction of such 
systems with the electromagnetic field may be treated in great detail to obtain valuable 
analytic results and, hopefully, the analysis of two-level systems generates insights that 
may be extended to more realistic situations. Fortunately, there are several important 
instances in which the application of the two-level model provides a very good 
approximation to a more complete theory. In the following, we label the upper level of the 
system by the letter a and the lower by the letter b. From Equation [ I- 12a ] we write, 
specifically, the wave function of the two level system as 

^r,t) = {r\\^f{t))= C^t) uj^r) exp{-i (Oj) + C,{t) u,(f) exp{-i(0,t) [ HI-l ] 

where we know from Equation [ I-12b ] that the time varying coefficients satisfy, in 
general, the following equations: 

Q(^)=-^{c,(?)(£,|i^i(?)|£,) + Q(?)<£,|i^,(OK>exp( [ni-2b] 

If we take the interaction to be the electric dipole interaction with an applied electric field we 
may write 

^i(0=~Ps"E(Rs>^)=-ef-E(R5,?) [in-3a] 
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where denotes the position of the center of the two-level system or atom.!"* Thus we 
(£, \!H it) I E\ =<U,.=-{E,\ei\E\-t(^,,t] [ ni-3b ] 



write 



In all but the most bizarre circumstances we may use persuasive symmetry arguments to 
reason that 

(£,|e?|£,>0 

Thus Equations [ 111-2 ] reduce to 

Cit) =- -^C^t) exp( mj) [ III-4a ] 

Clt) =- -jCXt) '^a^ exp(-/a),,?) [ ni-4b ] 

where V =-(£je x |^,)- e(r, =- p e(r, . 

rabi Flopping -- without damping; 

For an oscillatory applied field 

'^ah =-|pEjcosco^r=-^pE„exp(-ico^r) + c.c. [III-5 ] 

we see, in the rotating-wave approximation, that 



<^iO =^^E„C,(?) exp[ / (a),,-0)J?] 
= ^^0 C,(?)exp[/(a),,-a),)?] 



[ in-6a ] 



The use of this form of interaction needs considerable elaboration, but we defer that discussion until later. 
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An 

= \^fcXt) exp[-/(a),,-a),)?] 



[ ni-6b ] 



where Q^ = p E^/^ defines the so called Rabi flopping frequency which is, of 
course, a measure of the strength of the electromagnetic interaction. Clearly, the coupling 
terms have maximum effect when the frequency of the applied field is resonant with the 
level splitting. In most treatments the frequency detuning of the field is expressed as 
6a) = co^j -co^ and the system's wave function - i.e. Equation [ ni-1 ] - is written in terms 

of slightly modified time varying coefficients by transforming to the rotating frame of 
reference - viz. 




[ in-7a ] 




[ ni-7b ] 



,f) = CXt) exp[/ (-^Sco-CO J t\ ulr) + CXt) exp[/ (--^Sco-CO J t\ u^r) [ m-Vc ] 



The coupling terms in this rotating frame of reference become 




[ in-8a ] 




[ lll-8b ] 



or in matrix form 




[ III-8C ] 
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We look for a solution in the form C{t) = C{0) exp{^ I iCl^ tj where Q.^ is the 
generalization of the Rabi flopping frequency. Therefore, the condition 

det[M-^2'*I] = 



yields the generalized Rabi flopping frequency 

and the general non-damped time evolving wave function 



sm 



cos(in^) + ^sin{ir2^)j 



A(o) 



[ III-9a ] 



[ III-9b ] 



Rabi Flopping -- with Damping : 

Neglected interactions (e.g. spontaneous emission, collisions, thermal fluctuations) limit 
lifetime of a state of a two-level system. One class of lifetime limiting interactions may be 
described phenomenologically by adding decay terms to the equations of motion - i.e. 
to Equations [ III- 8 ] - as follows: 



or 



dt 



d 




( 
i ' 


dt 


(CM) 





-(5(0- /yJ 

ilf (5(0+ iY.) 



A(0 



[ m-lOa ] 

[ m-iob ] 



= ^M'C{t) [in-ioc] 
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Again we look for a solution in the fonn C[t) - C(0) exp( \ i At) so that 

det[M-Al] = 



or 



Therefore 



[in-lla] 



~ where the we refer to Yafe= ilYa"*" yJ average decay rate constant and 



n'*= J[6a)-it(Y«-Y.)]'+t^^| 



as the generalized complex Rabi flopping 

frequency — and the general time evolving wave function — in the rotation frame — may 
be written 



Ait) J 



cos 



(iU^t)- 



sm 



(lU^t) 



COS 



(lU-t) 



+ 



^si„(in«,) 

KY„-y„) + '8ra 



1 



sm 



(iU^t) 



exp(--2Yai t) 



'cXo) 

^Q(O) 



[m-iib] 



Density Matrix Treatment of a Two-Level Systems: 



Recall from Equation [ 11-30 ] 



n n 



Using |aX«l + I^X^I = 1 we may write 



P = -^{p(\a){a\ + \b){b\)^ -H (|a)(a|+|l;)(Z7|)p} [in-12a] 
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If 9{ = , then 

p=i{p[£ja)(a|+£,|&X^| + |aX«Ki + |^)(^K.] 

-[E^\a){a\ + E,\b){b\+fH:,\a){a\ + ^,\bXb\]p} 



Writing matrix elements ~ i.e. representatives ~ of the density operator 

{a\p\a)=-j,{{a\p\b){b\f^C,\a)-{a\i^,\b){b\p\a)} [IH-lSa] 

(Z7|p|Z7>=-^{(Z7|p|a)(a|itfJZ.>-<Z.|5/JaXa|p|^)} [ III- 13b ] 



(«|p|^>=-|{-(^«-^.)(«|p|^) + [(«|p|«)-(^|p|^)]<«Ki|^) 

+ {a\p\b)[(b\^ ,\b)-{a\^ ,\a)]} 



We may again assume, by symmetry arguments, that {a\ \a) = \b\ ^^\b)=0 so that the 
equations of motion for the elements of the density reduce to 

P..-P««= M"'2l^,,p,, + c.c.] [in-14a] 

Pab=-i(^abPab + i ab[p aa ' P bb] [ HI-Mb ] 

To introduce an element of reality, we add to these equations a pair of the most 
intuitively satisfying damping terms (a useful attribute of density matrix 
formulations) ~ viz. 

Pbb-ppaa=lifi~'^'^abPba+C-c]-^[(pbb-Paa)-(pbb-Paa)^] [HMSa] 
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[ ni-15b ] 



Transform these equations to a rotating frame by taking Pab~ Pab ^^P (~ ' ^) where 
is assumed to be a slowly varying function of time which satisfies the equations of 

motion 

p,, = -/6a)p,, + /rV^,exp( i(o^t) [p^-pj-yp,, [in-16a] 



Pbb- Paa = [i ^'^ab ^^r ^ P ba + ^c] 



-r 



{pbb-Paa)-(pbb-Paa)^] 



[ III-16b ] 



£xp (— ico^?) + c.c.j and ignoring terms proportional to exp(+ / 2(0 
-- i.e. the rotating wave approximation - we find that 

P,, = -/(5co-/Y)p,,-/(/2n^) [p„,-p,,] [in-17a] 

P - P «« = h' P + c-c] - r [(p - p - (p - p J [ iii-i7b ] 

Steady state behavior ; If we take all time derivatives in these equations equal to zero, 
we obtain 



Pat 



i^o [p - P _ ( 5co+ / y) [p - P 



[III-18a] 



r 



(Pbb-Paa)=(pbb-Paa). 



5a)'+Y' 



6a)'+Yll+f^„l7rY 



[ ni-18b ] 



or 
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-I 1 1 1 1 1 1 1 1 1 1 

2. 4. 6. S. 10. 

We have graphed Equation [ III- 18b ] at resonance - i.e. 5co= - and, from the following 
expression, we see how the oscillatory polarization is saturated at high electromagnetic 
powers. In general, the macroscopic polarization is then given byi^ 



That is to say 

\P> = (P P)aa +(P P).^, = ^aaPaa+ Kb Pba + Ka P ab + Kb Pbb = Kb Pba + Ka Pab 

since P^^ = P^^ = by symmetry. 
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If we graph Equation [ IE- 18b ] at resonance - i.e. 5a)= - we see how the oscillatory 
polarization is saturated at high electromagnetic powers - viz 

P = N (dipole moment of system) = A/^ p + c.c] [ III- 19a ] 

P= -T—^ — exp(jco, r) + c.c. [Ill- 19b] 

2(6a) +Y j 

N\^'[d(o-iy][p,,-pJ 
P= F T — ^E„exp(j(0, ?) + c.c, [111-19C] 

2^l5a)^+Y^[l+f2„t/rYj 
The Vector Model of the Density Matrix: 



There is a set of arguments by analogy which is exceedingly valuable in treating 
transient excitation problems in optics. The basis for the analogy lies in the fact that 
Equations [ 111-16] are identical in form to the famous Bloch equations of magnetic 
resonance. 1^ If we make a transformation to a rotating frame ~ i.e., 
M+ = M+ exp (+i co^ t) — the Bloch equations take on the form 



1^ Recall from the theory of magnetic resonance 

dM/dt= Y„,,MxH 

where Y^^„ is the gyromagnetic ratio. This equation of motion is greatly simplified if it is written in 
terms of the circular polarizations M+ = + i'M^ and H+ = + i — viz- 

Mi = +^-Y„.,M,H,±/Y™sM,H, 
and M^=i^ [M^H_ - M_Hj 

To these Bloch added the phenomenological longitudinal (thermal) relaxation time 7] and transverse 
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= +/ 5a) ± / Y„>ag exp (±/ (o^t)H^-U, /T^ [ in-20a ] 



M = 



H_ exp(-/ CO, t) + c.c. 



{M,-Mj/r, [ni-20b] 



where 5co =(0 , - y^.^^ . By comparing these Bloch equations with Equations [ III-16 ] 
and [ III- 17 ] we see that we have identical problems if we make the identifications 



i/T, = r 

l/T, = y 

In other words we have the equivalent equation of motion 

M= Mxh^ff- Y(M-zMj-rzM^ 

where the e^/ec^ive Rabi precession field ^2^^^ is given by 



^^ff= x-5ol) z 



[ III-21a ] 
[ ni-21b ] 
[ in-21c ] 
[ III-21d ] 
[ in-21e ] 
[III-21f ] 

[ III-22 ] 



[ ni-23 ] 



dephasing time T2 so that 



M, = +i y^ yuw+i y^^ H, - M,/r, 
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Initial state of system 



X 




'Rabi field" 



Effect of transient excitation 
- precession of state vector 




'Rabi field" 



State after 90° pulse 




"Rabi field" off 
-start of "free decay" period 
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